Twin Families of Bisolitons in Dispersion Managed Systems 
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Bisolitons in optical fiber lines with dispersion- 
management, were first discovered using computer mod- 
eling^ and later experimentally.^ Bisolitons can be 
viewed as a two component soliton molecule. In numer- 
ical simulations, they are stable over long propagation 
distances and, if perturbed, oscillate about equilibrium. 
We investigate the structure of such pulses assuming that 
fiber losses are completely compensated and propagation 
of pulses through optical fiber in a dispersion managed 
system is governed by the nonlinear Schrodingcr equa- 
tion: 

iuz + d{z)utt + ■y\u\'^u = (1) 

where u = u(t, z) is the slowly varying envelope of the 
electromagnetic field inside the fiber. We consider a sim- 
ple case of a piecewisc constant dispersion function d{z), 
where a fiber span of length Zdm/'^ with normal disper- 
sion alternates with equal-length spans of anomalous dis- 
persion fiber. The function d{z) can be represented as a 
sum of an oscillating part d{z) and a residual dispersion 
do such that d{z) = do+d{z). Here {d{z)) = 0; d{z) = di, 
if < z < Zdmf^, and d{z) = -di if Zd„/2 < z < Zdm- 
In this system, the characteristic length of the nonlin- 
earity is z„; ~ l/|Pp, where P is the peak power of the 
bisoliton, while the characteristic length of the residual 
dispersion is Zdo ~ T'^/do, where r is the pulse width. 
If the period of the dispersion map Zdm is much smaller 
than Zni and Zdo, then the spectrum u of the solution 
to Eq. is a slowly varying function of z on the scale 
Zdm <C Zni , Zdo and can be represented as 



q{oj, z) cxp —iuj'^ 



d{z')dz' 



(2) 



The exponential term captures the fast (in z) phase 
and q{uj, z) captures the slow amplitude dynamics of the 
spectral components. As has been shown, the evolution 
of the spectral components at leading order can be de- 
scribed by 



iqz{uj) - dauj'^q{uj) + -iR{q{uj),uj 
1 f sin(sA/2) 



0, where 

R{q{u)),uj) ^ J q{uJi)q{uJ2)q*{uJ3) 

xS{uJi + L02 — UJ3 — uj)dwiduj2duj3 (3) 



Here s = ZdmC?i/2 is dispersion map strength and A = 
Lol + — oj'l — Lo^ . We determine a shape of a bisoliton 
solution following earlier work by P.L.^ If a solitary wave 
solution with phase period has the form q{uj) = 
A[ijj)e^^^ , then amplitude A{uj) evolves according to the 
integral equation: 



-\A - doLO^A + 7i?(^(w), w) = 



(4) 



Rescaling variables t — tqt, uj = Q/tq and A{ll}) = 
aip{Q), where tq — s^^^, a = 27r{sX/^y^^, results in a 
dimensionless equation 

- ( 1 + don'^) ip{n) + R {(p{n), n) = O, where 
1 f sin(A/2) 



R{ipin),n) 



(27r)2 



A/2 



n)dnidn2dfi3, (5) 



which depends on a single parameter do = do/{sX). Here 

We study the structure of bisolitons as a function of 
do- Following the experimental work^ we consider an- 
tisymmetric solutions of Eq. ©. To solve this integral 
equation we use the iterative procedure: 



1 + doQ^ 



(6) 



Here Podd(/(a;)) = (/(x) - /(-x))/2 is a projection op- 
erator onto the set of odd functions and is an in- 
verse Fourier transform. A modified Petviashvili stabi- 
lizing factor^ Qn 



F-'[ip,,{n)] F- 



R{iPn{n),n) 



1 + dofl^ 



T=0.5 



allows the scheme to avoid trivial solutions ip = 0. The 
most costly part of this iterative procedure is evaluation 
of R, which involves a triple integral. To expedite evalu- 
ation of these integrals we used the procedure described 
by P.L.** It should be noted that the bisoliton solution 
of Eq.© represents the tmchirped pulse shape at the 
middle of each span with positive dispersion. 
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Fig. 1. Bisoliton energy as two valued function of dimen- 
sionless residual dispersion dg. Solid and dashed lines on 
insert correspond to the upper and lower branch bisoli- 
tons with do = 0.256. Horizontal and vertical axes on 
insert correspond to dimesionless time and amplitude. 



We begin by studying the solutions for the parame- 
ter do = 0.067 (a realistic value for communication sys- 
tems^), choosing (po as a sum of two shifted real valued 
Gaussian functions with opposite signs. The iteration 
procedure converges to the fixed point as the Petviashvili 
factor Qn approaches 1. The iteration is stopped when 
IQn — 1| < 10~^. The value of do was then varied in small 
increments. We use the solution found for the nearby do 
as the initial "guess" for the next value of do- Fig. ^ rep- 
resents bisoliton energy as function of do- Remarkably, 
this is a multiple valued function with two branches. Cal- 
culation of solutions on this second branch required solu- 
tion of the Arnoldi-Lanczos approximation problem for 
the linearization of iterration operator. Its limit point is 
located in the vicinity of c?oh/ — 0.426. With all other pa- 
rameters fixed, smaller values of do correspond to smaller 
values of residual dispersion. Therefore, the limit point 
corresponds to the largest value of do for which bisolitons 
are supported. According to our calculations, for values 
of do > dobf bisolitons will fail to exist and we will only 
observe a pair of interacting dispersion managed solitons 
that are not bound, not a bisoliton. The insert to the fig- 
ure shows that the higher energy bisoliton is wider, with 
greater separation and broader shape. 

Direct numerical simulations demonstrate stability of 
both the lower and the upper-branch bisoliton solutions 
over realistic distances (300 periods). For very long prop- 
agation distances, the upper branch showed signs of in- 
stability. Additionally, the pair of pulses composing the 
solution tends to stay bound whenever the pulses are 
pulled apart. The separation between the pulses spread 
apart in this manner oscillates about the separation for 
the bisolitonic solution. 

In the remainder of this paper, we will discuss the 
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Fig. 2. The contour plot of log(|(p|) as a function of time 
variable t and a parameter do- The dotted lines repre- 
sents the peak amplitudes of solutions. 



structure of the lower branch solutions. A later paper 
will present details about the upper branch solutions. 
The logarithmic density profile of lower branch solu- 
tions' amplitude for a range of values of do is shown 
in Fig. 121 There lighter shades of gray correspond to a 
higher value of the amplitude. The black lines correspond 
to zero values. The dashed lines indicate where the solu- 
tions have their maxima. A horizontal slice of this plot 
gives an amplitude profile for a fixed value of do- For 
example, for a value of do = 0.167 the logarithm of am- 
phtude is represented on Fig. El The dashed hue on the 
same graph represents the result of direct numerical sim- 
ulations of Eq. ^ after 300 dispersion map periods. A 
solution of equation jSl provided one boundary condition 
for this simulation (the launched pulse). 

As we see in Fig.|21bisoliton tails change sign for values 
of do < docr ( docr — 0.269 ). As the value of the di- 
mensionless residual dispersion becomes greater than the 
critical value do > doc i the phase of the tails remains 
unchanged, and the amplitude approaches an exponen- 
tially decaying function. Our equation Eq. Q reduces 
the many physical parameters from Eq. (Qto a single 
parameter. We have computed the bisolitons for a range 
of parameter values, and those computed solutions can 
be used to write solutions of the original physical system 




(7) 

where it was convenient to specify the solution at 
Zdm/^ + "tnzdm because it is chirp free at this point. 

We use u{t, Zdm/ ^) of Eq. Q as the initial condition 
in direct numerical simulation of Eq. ^ to study the 
dynamics of lower-branch bisoliton solutions for differ- 
ent values of do. In particular, we compare temporal- 
spatial behavior of solutions corresponding to values of 
docr < do < dob/ and < do < docr over a map period. 
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Fig. 3. Intensity profiles of a bisoliton with do — 0.167: 
solution of Eq. (jSJ (solid), and a result of propagation 
this solution over 300 periods through the line with s ~ 
1.5, do = 0.0125, and A = 0.05 (dashed) 



We consider a system with s = 1.5 and mean dispersion 
do = 0.0125. Fig. 4(a) and Fig. 4(b) represent propaga- 
tion of the initial pulse with do = 0.417 and do = 0.167 



respectively. For such choices of the system parameters, 
the values of phase periods A must be 0.05 and 0.02. 

The bisolitons for values of do above the critical value 
have no zeros other than at r = in their unchirped 
state, while if do < docr the solution will have an increas- 
ing number of zeros with smaller do . The example shown 



in Fig. 4(b) shows that for larger values of do the local 
minima which in part (a) are at z = split into pairs 
of minima, which are "shifted" towards Zdm/2. In fact, 
comparing the dynamics of do = 0.417 and dp = 0.167 
bisolitons indicates that the larger the value of dg the 
more the minima will shift from the narrowest states (the 
valleys) of the pulse to its broadest state (the ridges). 

In conclusion, we have used a slowly varying strobo- 
scopic equation to calculate bisoliton solutions with well 
resolved tails. This equation can be rescalcd so that it 
has a single dimensionless parameter dg. We have found 
a range of do such that there are two bisoliton solutions 
for each value of do. In addition, the structure of the tails 
for the lower branch solutions was described in terms of 
the value of do. 
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Fig. 4. Waveforms of a logarithm of the amplitude for 
bisoliton propagation through DM system with s = 1.5, 
do 0.0125, 7 = 1 (a) A = 0.05 and (b) A = 0.02. 
The soliton magnitude |m| is symmetric. Bisolitons are 
propagated over half of the period, from z^^/A to 3Z(i„/4 
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